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A  COMPARISON  OF  SMMR  AND  AIRBORNE 
MICROWAVE  MEASUREMENTS 
DURING  GOASEX 


1.0  INTRODUCTION 


The  Gulf  of  Alaska  Experiment  (GOASEX)  was  a  joint  NASA/NOAA 
effort  conducted  to  evaluate  the  performance  of  the  ocean 
surveillance  satellite,  SEASAT.  GOASEX  was  conducted  from  28 
August  to  26  September  1978,  shortly  after  the  SEASAT  launch  on 
26  June  1978.  Surface  and  near-surface  measurements  were  made 
by  a  variety  of  instruments  on  aircraft,  ships,  and  buoys  for 
comparison  with  measurements  made  by  the  remote  sensors  aboard 
the  satellite.  The  Naval  Research  Laboratory  RP-3A,  which 
carried  a  set  of  downward-looking  radiometers  to  measure  the 
microwave  emission  of  the  open  ocean,  participated  in  the 
experiment. 

The  SEASAT  remote  sensing  instruments  were  designed  to 
determine  various  oceanic  and  atmospheric  parameters  on  a  global 
all-weather  basis.  In  the  case  of  the  Scanning  Multichannel 
Microwave  Radiometer  (SMMR) ,  they  are  sea  surface  temperature, 
wind  speed,  and  atmospheric  liquid  water  and  water  vapor.  These 
parameters  are  determined  in  a  multi-step  process.  First,  the 
raw  SMMR  output  measurements  are  converted  to  the  absolute  power 
received  by  the  radiometer,  the  antenna  temperature,  TA.  Then 
T,,  is  corrected  for  the  directional  effects  and  antenna  losses 
to  obtain  the  absolute  radiation  from  the  scene  being  viewed  by 
the  radiometer,  the  brightness  temperature,  Tg.  Although  Tg  is 
a  weighted  average  over  the  instrument  parameters  of  bandwidth 
and  spatial  response,  it  is  otherwise  independent  of  the  sensor 
and  is  a  unique  measurement  of  the  scene  which  depends  only  on 
the  physical  properties  of  the  scene.  The  last  step  in  the 
process  is  the  conversion  of  Tg  to  the  desired  geophysical 
parameters  by  empirically-theoretically  determined  algorithms. 

Manuscript  approved  August  23,  1983. 
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The  most  obvious  method  to  evaluate  the  overall  instrument 
performance  is  to  see  how  accurately  it  determines  the  geophysi¬ 
cal  parameters  which  it  was  designed  to  measure.  Sea  surface 
temperature,  for  example,  can  be  measured  in-situ  and  compared 
to  the  corresponding  values  obtained  from  SMMR.  Instrument 
evaluation  by  this  method  involves  not  only  the  instrument 
itself,  but  also  the  algorithms  used  to  derive  the  parameters 
from  the  instrument  measurements.  A  more  direct  method  of 
validating  the  instrument  alone  is  to  validate  the  input  to  the 
algorithms,  i.e.,  the  brightness  temperature.  Of  the  quantities 
computed  during  the  data  analysis  process  of  converting  raw 
measurements  into  geophysical  parameters,  Tg  is  the  first  one 
which  is  independent  of  the  instrument  and  therefore  subject  to 
external  verification. 

This  verification  is  provided  for  SMMR  by  the  NRL  aircraft 
experiment:  radiometers  aboard  the  aircraft  measure  the  micro- 
wave  brightness  temperature  of  the  ocean  surface  for  comparison 
with  the  corresponding  values  measured  by  SMMR.  The  performance 
and  results  of  the  experiment  are  described  in  this  report. 

2.0  EXPERIMENT  DESCRIPTION 

2.1  AIRCRAFT 

Four  passive  microwave  radiometers  were  mounted  on  a  pallet 
in  the  lower  fairing  of  the  aircraft.  Figure  1  shows  the  RP-3A 
aircraft;  the  fairing  is  immediately  behind  the  front  landing 
gear.  One  of  the  radiometers  was  dual-frequency,  measuring  the 
horizontally  (H)  polarized  radiation  at  22.2  GHz  and  the  verti¬ 
cally  (V)  polarized  radiation  at  31.3  GHz.  The  other  three 
radiometers  were  single-frequency,  measuring  both  H  and  V  at 
6.15,  14.3,  and  19.3  GHz.  The  pallet  could  be  tilted  during 
flight  so  that  the  radiometers  could  view  either  in  the  nadir 
direction  or  at  50  degrees  off  nadir. 


The  NRL  research  aircraft  RP-3A  1^9670.  The  microwave 
radiometers  are  located  in  the  lower  fairing  behind  the  front 
landing  gear. 


All  of  the  radiometers  were  Dicke  switched.  The  6.15  and 
14.3  GHz  radiometers  used  tunnel  diode  amplifiers  covering  the 
bands  5.9  to  6.4  and  13.5  to  15.5  GHz.  The  other  three  were 
crystal-mixer-superheterodyne  receivers  with  local  oscillators 
at  19.3,  22.2,  and  31.4  GHz  and  IF  amplifiers  with  band  passes 
of  10  to  300  MHz.  The  22/31  GHz  radiometer  used  a  corrugated 
horn  with  beam  widths  of  8.4°  and  7.9°  and  main  beam  efficiencies 
of  97.5%  and  98.2%.  All  the  others  used  horn  lens  antennas  with 
beam  widths  and  beam  efficiencies  of  7.9°  and  80%,  7.9°  and  80%, 
and  8.0°  and  87%  at  6.15,  14.3,  and  19.3  GHz,  respectively. 

In  addition  to  the  radiometers,  the  aircraft  carried  a 
complement  of  instruments  to  make  supporting  meteorological  and 
environmental  measurements.  Wind  speed  and  direction  at  aircraft 
altitude  were  computed  by  the  aircraft  inertial  navigation 
system  (INS) .  Sea  surface  temperature  was  measured  by  a  Barnes 
infrared  radiometer.  Atmospheric  pressure,  outside  air  tempera¬ 
ture  (OAT) ,  and  dew  point  were  measured  at  aircraft  altitude. 

The  OAT  lapse  rate  and  water  vapor  scale  height  were  subsequently 
determined  from  these  data  taken  during  aircraft  descent. 

Whenever  possible,  environmental  parameters  measured  by  aircraft 
instruments  were  checked  against  independent  measurements  from 
ships  and  buoys. 

2.2  SMMR 

The  SEASAT  SMMR  instrument  is  an  imaging  microwave  radio¬ 
meter  measuring  both  V  and  H  polarized  radiation  at  the  five 
frequencies  6.6,  10.7,  18,  21,  and  37  GHz.  Half-power  beamwidth 
varies  from  0.9°  at  37  GHz  to  4.5°  at  6.6  GHz.  The  nadir  angle 
of  42°  at  the  794  km  satellite  altitude  results  in  an  incidence 
angle  of  48.8°  at  the  earth's  surface.  A  mechanically  scanning 
reflector  sweeps  the  line-of-sight  back  and  forth  across  a  swath 
of  600  km  width.  One  sweep  across  the  swath  takes  two  seconds, 
during  which  55  data  samples  are  taken  at  37  GHz,  14  samples  at 
6.6  GHz,  and  27  samples  at  each  of  the  other  three  frequencies. 
The  18  GHz  data  points  during  several  sweeps  are  shown  in 
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Figure  2,  illustrating  the  arc-shaped  sweep  and  the  region  of 
data  coverage.  Also  shown  are  the  resampled  grid  points  at 
which  the  Jet  Propulsion  Laboratory  (JPL)  determines  interpolated 
Tg's  (Njoku  et  al. ,  1980a).  A  detailed  description  of  the  SMMR 
instrument  and  its  method  of  operation  has  been  published  by 
Njoku  et  al.  (1980b) . 

2.3  FLIGHT  DESCRIPTION 

Five  GOASEX  flights  were  conducted  by  the  RP-3A  during  the 
period  19-27  September  1978.  There  were  equipment  problems  on 
the  flight  of  19  September,  but  good  data  were  obtained  on  the 
flights  of  21,  22,  25,  and  27  September.  This  report  is  based 
on  the  data  from  those  four  flights. 

The  corresponding  flight  paths  are  shown  in  Figures  3-6. 

The  first  two  flights  were  made  during  ascending  satellite 
passes,  and  the  second  two  during  descending  passes.  Each 
flight  began  with  a  transit  at  20,000  feet  altitude  from  the 
mainland  to  the  OSS  OCEANOGRAPHER.  Observations  began  during  a 
spiral  descent  over  the  OCEANOGRAPHER  to  a  500  ft.  altitude. 

Next,  a  star  pattern  was  flown  at  500  ft.  in  the  vicinity  of  the 
OCEANOGRAPHER.  The  RP-3A  then  flew  a  leg  out  to  the  subsatellite 
track,  and  a  leg  along  the  track  completing  the  measurements. 

The  radiometer  outputs  were  recorded  for  several  minutes 
alternately  at  0°  and  at  50°  incidence  angle  during  flight  by 
rotating  the  pallet,  so  that  several  minutes  passed  between 
successive  Tfi  determinations.  The  aircraft  altitude  varied  in 
steps  from  500  ft.  to  a  maximum  of  5000  ft.  during  the  two  legs, 
with  the  flight  generally  beginning  and  ending  at  500  ft.  The 
altitude  for  each  leg  is  shown  in  Figures  3-6.  The  total  measure¬ 
ment  period  lasted  approximately  3-4  hours.  Each  T_.  determination 

D 

corresponds  to  one  minute  or  more  of  continuous  observation, 
during  which  Tfi  generally  varied  by  less  than  2K.  Changes  in  Tg 
may  be  attributed  to  changes  in  environmental  conditions,  such 
as  cloud  cover,  surface  wind  speed,  or  aircraft  altitude. 
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LONGITUDE 


Figure  2. 


SMMR  data  swath,  illustrating  the  scans,  individual  18  GHz 
measurements  points,  and  JPL  resampled  grid  points. 


Figure  U.  NRL  RP-3A  flight  track  on  22  September  1978. 
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Figure  5*  NRL  RP-3A  flight  track  on  25  September  1978. 


OCEANOGRAPHER 
®  DATA  BUOY 

OCEAN  STATION  PAPA 


ALTIMETER  SUBSATELUTE  TRACK 
SMMR  SWATH  BOUNDARY 


CANADA 


% 

t  i; 


JMCC  JRD  AFB 

{"SITED  STATES 


3.0  DATA  ANALYSIS 


3.1  DATA  ANALYSIS  OUTLINE 

The  final  product  of  this  report  is  a  comparison  of  the 
brightness  temperatures  measured  by  SMMR  and  by  the  aircraft 
radiometers.  Since  the  two  sets  of  measurements  were  taken  at 
different  altitudes,  overlapping  but  slightly  different  frequen¬ 
cies,  and  slightly  different  nadir  angles,  several  calculations 
must  be  made  after  the  T  's  are  derived  before  a  direct  compari- 

Jj 

son  can  be  undertaken. 

There  are  four  main  steps  in  the  data  analysis  routine 
performed  by  NRL.  First,  brightness  temperatures  are  derived  at 
satellite  altitude  and  frequencies  from  SMMR  antenna  temperatures 
provided  to  NRL  by  JPL.  Second,  brightness  temperatures  at 
aircraft  altitude  and  frequencies  are  computed  from  the  flight 
measurements.  Third,  the  NRL  environmental  model  computer  program 
(Wisler  and  Hollinger,  1977)  is  used  to  compute  the  expected 
TB's,  atmospheric  loss  factors  and  emitted  atmospheric  radiation 
for  all  frequencies,  altitudes,  and  incidence  angles  for  the 
existing  environmental  conditions.  Finally,  estimates  for  the 
SMMR  T  '  s  are  derived  from  the  aircraft  T  's,  using  the  environ- 

D  D 

mental  model  calculations  to  correct  for  differences  in  frequencies, 
altitude,  and  incidence  angle  between  the  SMMR  and  aircraft 
measurements.  These  estimated  SMMR  T^'s  based  on  the  aircraft 
measurements  are  then  compared  with  the  SMMR  Tg ' s  measured  at 
SEASAT  altitude. 

The  four  data  analysis  steps  are  described  in  detail  in  the 
following  sections. 

3. 2  DERIVATION  OF  SMMR  BRIGHTNESS  TEMPERATURE 

SMMR  brightness  temperatures  at  satellite  altitude  are 
determined  from  antenna  temperatures  supplied  by  JPL  on  Sensor 
Data  Record  (SDR)  tapes.  The  SDR  tapes  contain  pairs  of  T^ 
values  for  all  frequencies,  and  earth  location  information  for 
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the  corresponding  arc-shaped  scans  (see  Figure  2) .  Computations 
made  at  NRL  pair  the  valu  j  for  H  and  V  polarization  and 
determine  the  corresponding  latitude,  longitude,  and  polariza¬ 
tion  angle,  a.  The  paired  TA  values  are  then  converted  to  T£  by 
JPL' s  interim  Antenna  Pattern  Correction  (APC)  method  rather 
than  by  their  final  detailed  and  extensive  APC  method  (Njoku  et 
al. ,  1980a) .  The  interim  method  assumes  a  relatively  homogeneous 
microwave  scene,  such  as  the  open  ocean,  and  uses  negligible 
computer  time  compared  to  the  final  method.  In  addition  to  the 
SDR  tapes  for  the  four  aircraft  flight  days,  JPL  has  provided 
NRL  another  tape  with  SMMR  TB  values  computed  at  JPL  by  both  the 
interim  APC  method  and  the  final  APC  method  for  one  flight  day, 

22  September.  Similar  data  were  not  provided  for  the  other 
three  flight  days  because  of  the  large  amounts  of  computer  time 
involved.  These  JPL  Tg's  are  used  to  check  the  accuracy  of  the 
Tfi  values  computed  at  NRL  from  the  JPL  SDR  tapes  by  the  JPL 
interim  APC  method.  First,  the  NRL  TB's  are  compared  to  the  JPL 
T  's  computed  by  the  interim  APC  method.  This  comparison  neces- 

D 

sitates  some  averaging  of  the  NRL  TB'S*  since  the  JPL  Tg's  are 
on  the  gridded  format  and  the  NRL  values  derived  from  the  SDR 
tapes  are  more  densely  packed  along  the  arc-shaped  SMMR  sweep 
(see  Figure  2).  Next,  the  NRL  TR ' s  are  compared  to  the  JPL  Tg's 
computed  by  the  final  extensive  APC  method,  and  a  small  bias 
between  the  two  Tn ' s  sets  is  determined  for  each  frequency  and 

D 

polarization.  The  bias  is  then  used  to  correct  the  NRL  TB 
values.  This  method  insures  that  the  NRL-computed  SMMR  TB's  are 
in  agreement  with  the  TB's  computed  at  JPL  by  the  most  accurate 
available  technique.  Also,  direct  conversion  of  the  SDR  tape 
T  ' s  to  T  ' s  results  in  a  greater  density  of  data  points  than  in 

H  D 

the  JPL  gridded  format,  and  thus  allows  the  selection  of  SMMR 
data  points  closer  to  the  aircraft  measurement  locations. 

A  brief  description  of  the  interim  APC  method  is  now 
given.  First,  the  antenna  temperatures  are  corrected  for  the 
cold  space  contribution  by 
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where  X  is  the  fractional  contribution  to  T  due  to  the  portion 
of  the  antenna  viewing  the  2.7K  space  background.  Values  of  X, 
provided  by  JPL,  are  listed  in  Table  1. 


Next,  Tb  is  derived  from  T^.  if  t^h  and  T^y  are  the  hori¬ 
zontal  and  vertical  antenna  temperatures,  and  and  T  „  the 

B<J>  B0 

horizontal  and  vertical  brightness  temperatures,  then  they  are 
related  by 


where 


TAH  GH0TB9  +  GH<rB<j> 


T  1  =  T  4.  r*  rp 

AV  v  0 1 B  0  vtj)  *B({) ' 


2  2 

GV0  =  gV9  COS  a  +  gV<p  sin  a  +  TERM  v 


Gv<p  =  gv<D  cos2a  +  gY0  sin2a  "  TERM  v 


3H0 


a  2  2 

^H0  cos  a  +  gH^  sin  a  +  TERM  H 


H<t> 


gH(j)  cos  a  +  gH0  sin  a  -  TERM  H 


TERM  V 

TERM  H 


4ve  9V(p  sin2a  cos  6 
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=  A 
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gH0  gH(J)  Sin  a  COS  6 
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The  above  equations  are  derived  from  the  method  of  Claassen  and 

Fung  (1974) ,  and  account  for  the  cross-polarization  correction 

due  to  the  instrumental  polarization  rotation  with  scan  ancle 
2  2 

(through  sin  a  and  cos  a),  cross-polarization  (through  the  g's) 
and  cross-track  gradients  (through  cos  6).  The  g's  are 
integrated  antenna  pattern  co-  and  cross-polarization  terms  and 
are  given  in  Table  1.  6  is  the  relative  phase  between  polariza¬ 

tion  components.  Values  of  cos  8  are  derived  empirically  at  JPL 
by  averaging  300  minutes  of  data  and  are  also  given  in  Table  1. 


SMMR  APC  Parameters 


37V  .015922  .9568  .0432 


Equations  (2)  are  inverted  to  obtain  Ton  and  T_,  from  T' 

DV  D(J)  H 

and  Tj^.  Note  that  since 

gve  +  gV4)  =  1  and  gH6  +  gHd>  =  1' 

then 


As  a  result,  the  inverted  equations  (2)  may  be  expressed  as 


gH4>tav 


—  P  ml 

V(j>  AH 


rn 

~B  4> 


p  ml  _  r>  ml 

v 0 1 AH  H6  AV 


(4) 


The  brightness  temperatures  computed  from  equation  (4) 
represent  the  JPL  interim  APC  method.  A  comparison  of  the 
resulting  Tfi  values  computed  at  NRL  from  the  SDR  tape  TA's  and 
the  interim  values  provided  by  JPL  shows  that  the  two  sets 
are  quite  close  for  all  frequencies.  The  average  RMS  difference 
is  0.7K,  showing  that  the  interim  APC  method  is  being  correctly 
applied  at  NRL.  On  the  other  hand,  the  final  JPL  Tfi  values  are 
a  few  degrees  lower  than  the  corresponding  interim  Tfi  values  on 
the  JPL-provided  tape  for  22  September.  The  difference  between 
the  Tg's  computed  by  the  final  APC  method  and  those  computed  by 
the  interim  APC  method  has  been  noted  by  JPL  (Christensen,  1981) . 
Therefore,  a  bias  is  derived  for  each  frequency  and  polarization 
between  the  NRL  Tg ' s  and  the  final  JPL  Tg ' s  by  averaging  both 
sets  over  the  entire  flight  area  for  22  September,  and  is  subse¬ 
quently  subtracted  from  the  individual  NRL  T_,  values.  The 
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biases  are  listed  in  Table  1.  Figures  7-11  show  the  resulting 
NRL  values  and  the  final  JPL  values  for  each  SMMR  frequency 
across  the  SMMR  swath.  Table  1  lists  the  RMS  deviation  between 
the  two  sets  for  each  channel.  The  low  RMS  values  of  approxi¬ 
mately  IK  indicate  that  the  SMMR  Tg  values  computed  at  NRL  with 
the  interim  APC  method  and  used  in  this  report  are  within  one 
degree  of  the  values  computed  by  JPL  using  the  final  APC  method. 


Figure  9*  A  comparison  of  lcS  GHz  SMMH  brightness  temperatures  derived  by 
JFL  and  by  NRL. 


es  derived  by 


The  bias  values  are  applied  to  all  four  flights,  which  occur  in 
the  same  area  during  a  six-day  period. 

3. 3  CALCULATION  OF  AIRCRAFT  BRIGHTNESS  TEMPERATURE 

Both  preflight  and  postflight  calibrations  of  all  aircraft 
radiometers  were  performed  in  the  laboratory  at  NRL  by  immersing 
the  antenna  in  black  body  enclosures.  The  preflight  calibrations 
occurred  within  one  week  of  the  first  flight,  and  the  postf light 
calibrations  within  one  week  of  the  last  flight.  Two  enclosures, 
spanning  the  range  of  measured  brightness  temperatures,  were 
used:  one  at  liquid  nitrogen  temperature  of  79K  for  the  low 
temperature  point,  and  the  other  at  ambient  temperature  of 
approximately  300K  for  the  high-temperature  point.  The  responses 
from  the  two  enclosures  are  used  to  linearly  interpolate  the 
equivalent  input  temperature  of  an  internal  noise  generator 
which  is  used  as  a  secondary  standard  for  inflight  calibration. 

Data  are  recorded  during  the  aircraft  flights  on  magnetic 
tape  and  on  strip  charts.  Antenna  readings  are  taken  alternately 
at  0°  and  50°  incidence  angle,  with  calibration  readings  being 
recorded  between  successive  pairs  of  antenna  readings.  One 
antenna  temperature  value,  T^,  is  calculated  for  each  antenna 
reading.  TA  is  then  converted  to  brightness  temperature,  Tg,  at 
aircraft  altitude  through  antenna  pattern  corrections.  The 
aircraft  APC  is  simpler  than  the  APC  for  SMMR,  because  the 
polarization  angle  a  is  always  zero,  and  also  because  beamfill¬ 
ing  is  more  complete  at  the  low  aircraft  altitude.  Conversion 
of  Ta  to  Tg  is  accomplished  by  inversion  of  the  equations 

TAH  =  g  TB0  +  {1_g)  TB0 

TAV  =  (1_g)  TR0  +  g  TB0, 

where  g  is  the  fraction  of  cross-polarized  power.  The  values  of 
g  are  computed  from  measured  antenna  patterns ,  and  are  approxi¬ 
mately  0.090  for  6.15  and  14.3  GHz,  0.055  for  19.3  GHz,  and 
negligibly  small  for  22.2  and  31.3  GHz. 
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The  aircraft  measurements  result  in  a  series  of  Tn  determin 

■D 

at ions  along  the  flight  path  and  within  the  SMMR  swath.  The 
first  Tg  location  is  the  star  pattern  area  in  the  vicinity  of 
the  OCEANOGRAPHER.  Several  of  these  measurements  at  50° 

D 

incidence  angle  are  selected  for  comparison  with  the  correspond¬ 
ing  SMMR  Tg  values. 

3.4  ENVIRONMENTAL  MODEL  CALCULATIONS 

Direct  comparison  of  the  SEASAT  SMMR  T  '  s  and  the  aircraft 
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Tg's  is  not  possible,  due  to  differences  in  frequency,  incidence 
angle,  and  altitude.  The  aircraft  measurements  must  be  correct¬ 
ed  to  correspond  to  the  SMMR  frequencies,  incidence  angle,  and 
altitude.  Differences  in  frequency,  angle,  and  altitude  are 
accounted  for  by  means  of  the  NRL  Environmental  Model  (Wisler 
and  Hollinger,  1977) .  The  model  is  used  to  calculate  the 
necessary  microwave  quantities  for  modifying  the  measured 
aircraft  Tg’s.  These  quantities  are:  (1)  the  radiation  emitted 
by  the  atmosphere  between  the  satellite  and  aircraft  altitudes, 

(2)  the  atmospheric  loss  factor  between  the  two  altitudes,  and 

(3)  the  brightness  temperature  at  aircraft  altitude. 

The  Environmental  Model  calculations  require  the  input  of 
values  for  several  atmospheric  and  oceanic  parameters  which 
largely  determine  the  microwave  brightness  temperature,  such  as 
ocean  surface  temperature,  atmospheric  water  vapor  and  liquid 
water  content,  low  altitude  wind  speed,  etc.  These  values  are 
derived  from  measurements  made  during  the  experiment  by  aircraft 
ships,  and  satellites.  Table  2  lists  the  input  parameters  and 
their  measurement  sources.  The  OCEANOGRAPHER  measurements  are 
published  in  the  GOASEX  Surface  Truth  Data  Inventory. 
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TABLE  2 


Sources  of  Environmental  Data 


Environmental 

_ Parameter _ 

Surface  Pressure 

Air  Temp,  at  Surface 

Temp.  Lapse  Rate 

Absolute  Humidity  at  Surface 

Water  Vapor  Scale  Height 

Clouds 

Near-Surface  Wind  Speed 
Sea  Surface  Temperature 


_ Source _ 

OCEANOGRAPHER 
OCEANOGRAPHER 

Aircraft  OAT  During  Descent 
OCEANOGRAPHER 

Aircraft  Data  During  Descent, 
Wentz  Data 

OCEANOGRAPHER,  Aircraft 
Observations,  Satellite  Imagery 

OCEANOGRAPHER,  Aircraft  INS 

OCEANOGRAPHER,  Aircraft  PRT-5 


For  each  location  along  the  aircraft  flight  path  at  which  a 
SMMR/aircraft  Tn  comparison  is  to  be  made,  environmental  model 
calculations  are  performed  using  the  best  set  of  environmental 
parameter  values  for  that  particular  location  and  time.  Some  of 
the  parameters  are  measured  or  determined  only  once  on  each 
flight  and  are  assumed  constant  for  that  flight.  For  example, 
the  temperature  lapse  rate  is  determined  from  OAT  measurements 
during  aircraft  descent  from  20,000  ft.  to  500  ft.  altitude 
immediately  prior  to  flying  the  star  pattern. 


The  most  difficult  parameter  to  determine  is  cloud  cover, 
i.e.,  atmospheric  liquid  water  density  versus  altitude  and  earth 
location.  To  assist  in  this  problem,  GOES  satellite  IR  images 
are  obtained  from  the  NOAA/Satellite  Data  Services  Division,  and 
the  images  are  interpreted  by  a  NOAA  meteorologist  to  provide 
estimates  of  cloud  type  and  thickness  and  of  percent  cloud 
cover.  Figures  12-15  show  the  IR  images  for  the  four  flights, 
together  with  the  corresponding  flight  paths.  Cloud  cover  in 
the  flight  area  is  generally  heavier  on  the  second  two  flights 
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Figure  13.  GOES  satellite  infrared  image  of  clouds  in  the  GOASEX  area  on 
22  September  1978  with  the  aircraft  flight  track. 


►  I 
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than  on  the  first  two.  The  information  from  the  IR  images  is 
used  together  with  cloud  information  from  the  other  sources  to 
estimate  the  atmospheric  liquid  water  content. 

To  obtain  a  valid  brightness  temperature  comparison,  data 
should  be  taken  only  from  regions  in  which  the  meteorological 
conditions  result  in  a  relatively  homogeneous  microwave  picture, 
free  of  strong  local  T_  gradients.  Otherwise,  the  actual  T0 
within  the  small  spot  size  of  the  low-altitude  aircraft  radiom¬ 
eters  could  be  sianif icantly  different  from  the  average  T_, 
within  the  large  spot  size  of  the  high-altitude  SMMR  radiometers. 
For  the  purpose  of  checking  on  the  degree  of  T_  variability, 
microwave  images  of  the  GOASEX  region  are  derived  from  the  SMMR 
37  GHz  H  measurements  and  are  shown  in  Figures  16-19.  The 
location  of  the  flight  path  in  each  image  indicates  that  the 
aircraft  generally  flew  through  quiet  microwave  regions  with  T0 
<  165K.  Regions  of  strong  microwave  variability  and  T_  >  180K 
correspond  to  heavy  cloud  areas  seen  in  the  IR  images.  One 
exception  is  that  the  areas  of  high  in  the  top  of  Figure  18 
represent  the  Canadian  landmass  rather  than  clouds. 

3.5  DERIVATION  OF  ESTIMATED  SMMR  BRIGHTNESS  TEMPERATURE 

Estimates  of  the  SMMR  Tg ' s ,  based  on  the  measured  aircraft 
Tg's,  are  derived  for  comparison  with  the  measured  SMMR  TB ' s 
according  to  the  following  scheme. 


A.  Look  Angle  -  The  measured  aircraft  Tg ' s  are  extrapo¬ 
lated  in  look  angle  from  the  50°  aircraft  angle  to  the  48.8° 


SMMR  angle.  An  estimate  E(FA,  48.8,  low)  of  Tg  at  aircraft 
frequency  Fa,  48.8°  incidence  angle  and  aircraft  altitude  is 
derived  from  the  corresponding  measured  value  M(FA,  50,  low)  by 
calculating  the  difference  in  T_  due  to  the  angular  change  and 
adding  the  difference  to  the  measured  value: 


E (Fa ,48.8, low)  =  M(FA,50,low)  +  [C (Fft , 48 . 8 , low) 

-  C(Fa,50,1ow) ) 


(5) 
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Figure  18.  Contour  map  of  SMMR  T  at  37  GHz  H  polarization  in 
the  GOASEX  area  on  25  September  1978  with  the 
aircraft  flight  track. 
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Figure  19.  Contour  map  of  SMMR  TB  at  37  GHz  H  polarization  in 
the  GOASEX  area  on  27BSeptember  1978  with  the 
aircraft  flight  track. 
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where  C  refers  to  the  value  of  TR  computed  from  the  NRL  Environ¬ 
mental  Model. 


The  calculated  difference  in  T_  is  2  K  or  less. 

fcs 


B.  Frequency  -  The  results  of  (A)  are  used  to  obtain  an 
estimate  E  (Fc  ,  48 . 8  ,  low)  of  T_,  at  SMMR  frequency,  F„  ,  and  air-- 

b  £5  b 

craft  altitude.  Figure  20  shows  an  example  of  computed  T 
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versus  frequency  for  both  V  and  H  polarization  with  the  location 


of  all  SMMR  and  aircraft  frequencies.  The  T_  curves  are  computed 
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for  aircraft  altitude,  48.8°  incidence  angle,  and  for  meteorolog¬ 
ical  conditions  of  the  22  September  flight,  and  are  shown  simply 
to  illustrate  the  relative  SMMR  and  aircraft  frequencies  and  the 


expected  T  differences  between  the  two  sets  of  frequencies. 
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Obviously,  T0  at  a  given  F0  can  be  derived  sometimes  by  extrapo- 
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lation  from  a  nearby  F^  and  sometimes  by  interpolation  between 


two  frequencies  F^  and  F^  •  For  extrapolation. 


E (Fg ,48.8, low)  =  E (F^, 48.8, low)  +  [C (Fg , 48 . 8 , low) 


(6A) 


-  C(FA,48.8,low)] 


For  interpolation. 


E(Fs,48.8,low)  =  fECF^^.S^Ow)  +  (1-f)  E  (Fft2 , 4  8 . 8  ,  low)  (6B) 

where 


f  = 


C (Fa2 ' 48 " 8 ' low^  ”  c (Fs • 48 • 8 f low) 


C(FA2'48‘8,1ow)  “  C(Fa1,48.8,1ow) 


c*  Altitude  -  The  results  of  (B)  are  extrapolated  to  the 
satellite  altitude,  i.e.,  the  low  altitude  estimate 
E (Fc ,48.8, low)  is  used  to  derive  the  high  altitude  estimate 
E(F  ,48. 8, high)  for  direct  comparison  with  the  SMMR  measurements: 

E (high)  =  L ( low, high)  E(low)  +  T  (low,  high)  (7) 

where  L  is  the  atmospheric  loss  factor  and  Tup  is  the  upwelling 
atmospheric  radiation  between  the  two  altitudes.  The  arguments 
F  and  48.8,  which  are  common  to  all  terms  in  equation  (7),  have 
been  dropped.  L  and  T  are  calculated  by  the  Environmental 


FREQUENCY,  GHz 

Figure  20.  Sample  brightness  temperature  curves  illustrating 
the  relative  position  of  the  SMMR  frequencies  and 
the  aircraft  frequencies. 


Model.  The  resulting  estimated  increase  in  T_  due  to  altitude 
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extrapolation,  E  (high)  -  E  (low) ,  varies  from  approximately  2  K 
at  6.63  GHz  to  12  K  (clear  sky)  or  21  K  (cloudy  sky)  at  37  GHz. 

When  equations  (5)  -  (7)  are  combined,  we  obtain  the  esti¬ 
mate  E (Fg,48. 8, high)  as  a  linear  function  of  the  measurements 
M (FA, 50 , low) .  Only  one  measurement  is  involved  when  frequency 
extrapolation  is  used, 

E(Fgj ,48. 8, high)  =  a1M (FA, 50 , low)  +  a2,  (8 

but  two  measurements  are  necessary  for  interpolation, 

E(Fgj ,48. 8, high)  =  a^M (FA1 , 50 , low)  +  b2M (FA2 , 50 , low)  +  b.^.  (9 

The  terms  a^,  a2 ,  b^,  b2 ,  and  b^  depend  only  on  the  Environmen¬ 
tal  Model  calculations  and  are  totally  independent  of  the 
measured  aircraft  brightness  temperatures. 

As  seen  in  Figure  20,  most  of  the  SMMR  frquencies  are  close 
enough  to  one  of  the  aircraft  frequencies  that  the  Tg  estimate 
can  be  extrapolated  from  a  single  aircraft  measurement  according 
to  equation  (8).  An  exception  is  the  SMMR  frequency  10.69  GHz, 
which  lies  intermediate  between  the  aircraft  frequencies  6.15 
and  14.3  GHz.  Consequently,  the  Tg  estimate  for  10.69  GHz  is 
derived  from  the  aircraft  measurements  at  6.15  and  14.3  GHz  by 
interpolation  according  to  equation  (9) .  Also,  the  SMMR  37  GHz 
is  far  removed  from  the  closest  aircraft  frequencies  of  22.2  GHz 
(H  polarization)  and  31.3  GHz  (V  polarization),  so  that  there 
will  be  more  uncertainty  in  the  resulting  Tg  estimates  at  37  GHz 
than  those  at  any  of  the  other  SMMR  frquencies.  Table  3  lists 
the  SMMR  frequencies,  the  corresponding  aircraft  frequencies, 
and  the  methods  used  in  deriving  the  T_  estimates. 


TABLE  3 


SMMR  and  Aircraft  Frequencies  and  Polarization 


SMMR 


6.63  v 
6.63  H 

10.69  V 
10.69  H 

18  V 
18  H 

21  V 
21  H 

37  V 
37  H 


Extrapolation 


Interpolation 


Extrapolation 


Extrapolation 


Extrapolation 


Aircraft 


6.15  V 
6.15  H 

6.15  V  and  14.3  V 
6.15  H  and  14.3  H 

19.3  V 
19.3  H 

19.3  V 

19.3  H 

31.3  V 
22.2  H 


4.0  RESULTS 

A  total  of  eight  separate  locations  are  selected  from  the 
four  flights  for  T0  comparisons.  These  locations  include  the 
four  star-pattern  areas.  The  other  four  locations  are  chosen 
near  the  center  of  the  SMMR  swath  where  polarization  angle  a  is 
close  to  zero  and  cross-polarization  errors  are  consequently 
minimized.  Table  4  lists  for  each  location  the  latitude, 
longitude,  date,  and  time  at  which  the  aircraft  measurements 
were  made.  For  the  star  patterns,  which  took  approximately  one 
hour,  the  mid-time  and  mid-location  are  given. 

TABLE  4 

SMMR/Aircraf t  T_  Comparison  Locations 

£5 


Date 


21  Sept 

1 

8 

21  Sept 

1 

8 

22  Sept 

1 

6 

22  Sept 

1 

7 

22  Sept 

1 

8 

25  Sept 

0 

6 

27  Sept 

0 

6 

27  Sept 

0 

8 

North 
Lat . 


East 
Long . 


00 

49.0 

218.3 

54 

48.4 

217.2 

50 

48.7 

226.3 

36 

48.9 

225.3 

27 

49.7 

222.3 

20 

48.9 

226.8 

50 

50.3 

215.3 

10 

48.5 

218.7 

SEASAT 
Pass  No. 


SEASAT 

Overpass 

Time 


The  brightness  temperatures  at  aircraft  altitude  are 
averaged  for  the  eight  locations  and  plotted  versus  frequency  in 
Figures  21-23  for  H  polarization,  V  polarization,  and  the  sum  of 
H  and  V.  Table  5  lists  the  averaged  values  for  H  and  V.  Theore 
tical  brightness  temperatures  are  calculated  for  each  case  by 
the  Environmental  Model,  averaged,  and  also  shown.  Two  sets  of 
environmental  calculations  are  made,  one  set  assuming  no  clouds 
and  the  other  set  assuming  complete  cloud  cover,  i.e.,  represent 
ing  minimum  or  maximum  cloud  effect  on  Tn .  Figures  21-23  show 
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that  there  is  reasonable  agreement  between  the  aircraft 
measurements  and  theory.  An  exception  is  that  there  is  better 
agreement  for  H  +  V  than  for  H  or  V  individually  at  6.15  and 
14.3  GHz,  indicating  the  possibility  of  some  depolarization 
error  at  these  two  frequencies. 

TABLE  5 

Brightness  Temperatures  at  Aircraft  Altitude 


Aircraft  TB,  K 


Frequency,  GHz 

H 

V 

6.15 

88.2 

148.8 

14 . 3 

98.9 

156.9 

Is.  3 

111.4 

175.6 

22.2 

141.7 

- 

31.4 

- 

194.1 

The  eight  sets  of  aircraft  measurements  are  then  extrapo¬ 
lated  to  SMMR  angle,  frequency,  and  altitude  according  to  the 
previously  outlined  method  to  obtain  esti>  ates  of  SMMR  Tg 1 s . 

The  calculations  are  performed  twice,  assuming  first  no  clouds 
and  then  assuming  clouds.  This  method  results  in  a  spread  in 
the  estimated  SMMR  Tg's.  The  actual  SMMR  measurements  should 
fall  within  this  range,  since  the  large  SMMR  footprints  are 
likely  to  be  partly  cloud-filled  and  partly  clear.  Figures  24- 
26  show  the  averaged  estimated  SMMR  Tg's  for  H  polarization,  V 
polarization,  and  the  sum  of  H  and  V.  Note  that  the  relation 
between  the  aircraft  measurements  and  theory  at  aircraft  altitude 
seen  in  Figures  21-23  is  preserved  in  the  relation  between  the 
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22.  Measured  aircraft  Tg  and  calculated  Tg  for  V 

polarization  at  aircraft  altitude,  averaged  over 
eight  locations. 
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Figure  23.  Measured  aircraft  Tg  and  calculated  TB  for  sum  of  V 
and  H  polarizations  at  aircraft  altitude,  averaged 
over  eight  locations. 


estimated  SMMR  Tg's  and  theory  at  satellite  altitudes  seen  in 
Figures  24-26.  This  result  indicates  that  the  extrapolation  of 
the  aircraft  measurements  from  aircraft  conditions  to  satellite 
conditions  is  being  performed  correctly. 

Next,  the  actual  SMMR  measurements  for  the  eight  locations 
are  determined,  averaged,  and  plotted  in  Figures  24-26.  The 
SMMR  Tg's  in  the  star  areas  are  determined  by  averaging  all  SMMR 
measurements  occuring  within  the  region  covered  by  the  RP-3A  in 
flying  the  star  pattern.  For  the  non-star  locations,  the  SMMR 
measurements  are  averaged  in  square  areas  around  the  locations, 
the  size  of  the  square  being  20  km  for  37  GHz,  60  km  for  6.63 
GHz,  and  35  km  for  all  other  frequencies.  The  sizes  of  the 
squares  are  chosen  so  that  a  minimum  of  two  SMMR  measurements 
are  averaged  at  each  location  for  each  frequency. 

Figures  24-26  show  that  the  SMMR  Tg's  are  generally  low 
compared  to  the  estimates  derived  from  the  aircraft  measure¬ 
ments,  except  at  6.63  GHz  H  and  V  and  at  21  GHz  H.  The  aircraft 
measurements  and  the  theory  are  comparatively  consistent, 
especially  in  the  H  +  V  curve  of  Figure  26,  where  depolarization 
errors  are  nullified.  It  must  be  concluded  that  there  are 
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significant  errors  in  the  measured  SMMR  brightness  temperatures 
at  most  frequencies. 

Table  6  summarizes  the  data  shown  in  Figures  24  and  25;  the 
aircraft  TD  is  the  average  of  the  two  values  for  clouds  and  for 
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no  clouds.  The  resulting  bias,  listed  in  the  last  column  of 
Table  6,  is  the  amount  which  must  be  added  to  the  SMMR  Tg  to 
bring  it  into  agreement  with  the  aircraft  Tfi .  Most  of  the  bias 
values  are  greater  than  5  K,  i.e.,  signifies.. .  ly  large. 
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Figure  25.  Measured  SMMR  T„,  estimated  SMMR  from  aircraft 
measurements,  and  calculated  TB  for  V  polarization 
at  satellite  altitude,  averaged  over  eight  locations 
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Figure  26. 


Measured  SMMR  T„,  estimated  SMMR  T_  from  aircraft 
measurements,  and  calculated  TR  for  sum  of  V  and  H 
polarization  at  satellite  altitude,  averaged  over 
eight  locations. 


TABLE  6 


Brightness  Temperatures  from  SMMR  and  Aircraft  at  SFASAT  Altitude 

Extrapolated 

SMMR  Tg,  K  Aircraft  Tg,  K  Bias,  K 


6.63  H 

92.0 

92.5 

0.5 

6.63  V 

149.2 

148.7 

-0.5 

10.69  H 

88.4 

98.2 

9.8 

10.69  V 

149.4 

152.8 

3.4 

18  H 

107.4 

114.7 

7.3 

18  V 

160.2 

173.8 

13.6 

21  H 

143.0 

143.8 

0.8 

21  V 

184.6 

193.7 

9.1 

37  H 

145.9 

160.6 

14.7 

37  V 

194.8 

2C9.2 

14.4 

The  SMMR  Tg's  are  used  as  inputs  to  various  geophysical 
algorithms  for  the  retrieval  of  geophysical  parameters  such  as 
sea  surfaces  temperature  (SST)  and  ocean  wind  speed.  Incorrect 
Tg  values  will  result  in  the  production  of  incorrect  geophysical 
parameter  values  from  the  algorithms.  Two  methods  have  been 
used  to  compensate  for  this  problem:  (1)  adding  a  predetermined 
bias  to  correct  the  SMMR  Tg's  before  using  them  in  the  algorithm, 
or  (2)  adding  a  predetermined  bias  to  correct  the  geophysical 
parameter  after  retrieval  by  the  algorithm.  In  the  Wentz  geophy¬ 
sical  algorithm,  the  first  method  is  used.  The  SMMR  TD  biases 

D 

derived  here  and  listed  in  Table  6  may  be  compared  to  recent 
values  used  in  the  Wentz  algorithm  (Wentz,  F.J.,  private  communi¬ 
cation).  The  two  sets  of  biases  are  plotted  in  Figure  27.  They 
are  in  reasonable  agreement  and  the  frequency  dependence  is 
identical.  An  exception  is  that  the  NRL  bias  values  at  37  GHz 
are  much  greater  than  the  corresponding  Wentz  values,  possibly 
because  the  NRL  values  are  based  on  measurements  at  the  much 
lower  frequencies  of  22  and  31  GHz. 

The  results  of  this  report  provide  independent  experimental 
verification,  both  qualitatively  and  quantitatively,  of  the 
empirical  observation  that  the  brightness  temperatures  measured 
by  the  SEASAT  SMMR  are  in  error  at  most  frequencies. 
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